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Abstract: Generating functions encoding cubic interactions of (partially-)massless higher- 
spin fields are provided within the ambient-space formalism. They satisfy a system of 
higher-order partial differential equations that can be explicitly solved due to their factor- 
ized form. We find that the number of consistent couplings increases whenever the squares 
of the field masses take some integer values (in units of the cosmological constant) and 
satisfy certain conditions. Moreover, it is shown that only the supplemental solutions can 
give rise to non-Abelian deformations of the gauge symmetries. The presence of these 
conditions on the masses is a distinctive feature of (A)dS interactions that has in general 
no direct counterpart in flat space. 
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1 Introduction 



Since the gravitational-interaction problem of higher-spin (HS) fields 1 was overcome by 
turning on a cosmological constant [8, 9], many studies have been devoted to understand 
HS field theories in (A)dS background. The most important result of these efforts are 
Vasiliev's equations [10, 11], that, together with String Theory, represent the only known 
frameworks in which interactions of HS particles can be consistently described. However, 
a deeper understanding of them requires further investigations, especially in relation to 
the possibility that they be only particular members of a wider class of consistent HS 
theories. A preliminary step towards this goal would be to construct the most general 
consistent cubic couplings. 2 Among them, only a subset is expected to be compatible with 
higher-order interactions [31-33] and lead eventually to fully non-linear theories. 

Finding HS interactions in (A)dS has been the aim of our previous investigations 
[34-37], where, making use of the ambient-space formalism [38, 39], we provided the trans- 
verse and traceless (TT) parts of all possible cubic interactions involving massive and 
massless totally-symmetric HS fields. The key point of our construction was to recast 
the consistency conditions for the cubic vertices into a system of linear partial differen- 
tial equations (PDEs), 3 whose solutions are in one-to-one correspondence to the consistent 
couplings. However, while for the case of massive and massless fields one gets relatively 
simple second-order PDEs, when partially-massless (PM) fields 4 enter the interactions, due 
to their higher-derivative gauge transformations, one has to solve higher-order PDEs. In 
[35] we only provided some examples of this kind of interactions relying on a numerical 
algorithm. 

This paper is aimed at completing the program initiated in [35] and at deriving the 
generating functions of cubic interactions involving PM fields. A key point of our result is 
that the number of consistent couplings depends, in a non-trivial way, on the masses of the 
interacting fields. More precisely, whenever the squares of the masses take some integer 
values and satisfy certain conditions, a new class of solutions appears. Since the solutions — 
which are not of the latter type — do not lead to deformations of the gauge symmetries, the 
aforementioned conditions on the masses represent necessary conditions for the presence 
of non-Abelian interactions. After some preliminaries on the ambient-space formulation of 
HS interactions, we summarize our results in Section 1.2. 

1 See e.g. [1, 2] for recent reviews on HS field theories. See [3, 4] for some reviews on Vasiliev's equations, 
and [5-7] for AdS/CFT-related issues. 

2 Many efforts have been devoted in this direction: see the references in [1] for an exhaustive list of 
works, and in particular the latest works [12-23] and [24-30] for flat-space and (A)dS interactions involving 
totally-symmetric HS fields, respectively. 

3 See the partial result [40] for the generalization of this method to mixed-symmetry HS interactions. 

4 The PM spectrum has been first discovered for lower-spin (spin 2 and 3/2) fields in [41-47], while its 
HS generalization has been considered in [48-51]. The implications of PM fields to (A)dS/CFT have been 
discussed in [52, 53]. In [54-59], various formulations for the description of PM fields have been proposed. 
See [60-63] for the interactions of PM fields, and [63-67] for their connection to conformal theories and 
massive gravity. Finally, see [68] for the representations of 50(1,4) . 
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1.1 Preliminaries 



Ambient-space formalism A way of describing totally-symmetric (A)dS tensor fields, 
'Pfti—Ha ! is through ambient-space fields, <& Ml --M a > that are subject to the homogeneity and 
tangentiality (HT) conditions: 

Homogeneity : (X ■ d x - U ■ d v + 2 + p) §(X, U) = , 
Tangentiality : X ■ d v U) = , (1.1) 

where &(X, U) is the generating function of <& Mi ---m s '■ 

oo 1 

*(X, U) = J2~ $m,..m s (X) U M - ■ ■ ■ U M ° . (1.2) 
When [i = 0, l,...,s — 1, the constraints (1.1) allow higher-derivative gauge symmetries: 

$(x, u) = (u- d x y +1 n(x, u) , (i.3) 

with gauge parameters Q satisfying 

(X -d x -U -du - n)Q(X,U) =0, X-duQ(X,U) = 0. (1.4) 

Massless fields, /u = , are the first members of a class of short representations where 
the other members, p = 1, 2, . . . , s — 1 , are PM fields. For other values of yu, no gauge 
symmetry is allowed, implying that the corresponding fields are massive. Focussing on 
unitary representations, the corresponding values of (i are constrained to the regions shown 
in Figure 1. Let us mention that, contrary to massless fields, PM representations are unitary 
only in dS. 

Im(/i) 



AdS 


dS 


( 


) 1 


2 


s-2 s-1 


8-1 + ^2 



Re(/i) 



Figure 1. Unitary values of /j, (for s > 0) . 



Cubic-interaction problem The most general expression for the TT parts of the cubic 
vertices reads 

S (3) = \. I C(Y, Z) $(X l5 Ui) $(X 2 , U 2 ) $(X 3 , U 3 ) Xi=x , (1.5) 

6 - J(A)dS Lr i= 

where /^)ds * s an integral over the codimension-one hypersurface X 2 = e L 2 (e = 1 for dS 
and e = — 1 for AdS), while C is an arbitrary function of the following parity-preserving 
Lorentz invariants: 5 

Yi = d Ui ■ d Xi+1 , Zi = d Ui+1 - du^ [i ~ i + 3] . (1.6) 



5 Henceforth, for brevity, we shall denote the dependence on all six variables by Y and Z- 
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Assuming the i-th field to be (partially-)massless ((P)M) (i.e. Hi € N), the correspond- 
ing compatibility condition of the cubic vertices (1.5) with the gauge symmetries (1.3) is 
equivalent to imposing 

[C(Y,Z),(U i -d Xi r +1 ] rr _ n ^0. (1-7) 
Using Leibniz's rule, one can recast the condition (1.7) into a higher-order PDE: 

Ci(fii - Hi) Ci(fii - m + 2) ■ ■ ■ £i(fii + Hi) C(Y,Z) = [pi := Hi-i ~ fH+l] , (1-8) 
consisting in the product of the following commuting differential operators: 

Ct(x) := Y i+1 d Zi _ x - Y^ x d Zi+1 + | (y i+1 d Yi+1 - y,_! dy^ + \ x) d Yi . (1.9) 

Depending on the number of (P)M fields involved in the interactions, one can have up to 
three PDEs, whose solutions encode all possible consistent couplings. 

1.2 Summary of the results 

As already mentioned, a peculiar feature underlying (A)dS interactions is the appear- 
ance of non-trivial conditions on the mass values for which the number of consistent cou- 
plings may get enhanced. More precisely, if the i-th field is at one of its (P)M points 
Hi £ {0, 1, . . . , Si — 1} , the solution space of the corresponding system of PDEs may be- 
come bigger whenever the conditions 

Hi + Hi+i ~ l-H-i G 2Z, (1.10) 

hold. An explicit analysis shows that, while for arbitrary values of the Hi' s the solutions 
can be written in terms of arbitrary functions of the following operators: 

Hi := d Xi+1 ■ dx^ d Ui+1 ■ du^ - du^ ■ d Xl+1 d Ui+1 ■ dx^ , (1.11) 

when the condition (1.10) is satisfied, additional solutions involving the operator 

G{Y, Z) := Y 1 Z x + Y 2 Z 2 + Y 3 Z 3 , (1.12) 

also appear. However, some of the G-couplings can be also expressed in terms of the 
-ff-couplings, so that the two kinds of solutions may have some overlap. Let us stress 
that, since the i/-couplings are trivially gauge invariant, non-Abelian interactions are only 
among those G-couplings which cannot be written as i7-couplings. Therefore, a necessary 
condition for non-Abelian interactions to be present is (1.10). The existence of the lat- 
ter mass-pattern may have also some interesting consequences for the interactions of one 
massless and two massive fields: on the mass-pattern, indeed, such interactions can induce 
deformations of the gauge symmetries related to non-trivial Noether currents involving 
fields with different masses, a novelty of (A)dS interactions. 6 

Our results are summarized in the following framed paragraph. 



6 This point has been omitted in our previous work [35], as, in the analysis of the flat limit (see Appendix 
D), we overlooked the singular points of the PDEs. 
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(P)M Massive Massive 



For arbitrary /J>2 — H3 : 



Ml 



^ = E ^ ^ ( y2 ' Zi,H 2 ,H 3 ) 

cti=0 

For //1 + yu,2 — Hz £ 2 Z , one also has 

c= £ z; 1 z; 2 z 3 *» y 2 

(-ri,T2,T3)e^fi 

x e~i v K T1T2TS (Y,G) 

G=G(Y,Z) 

where R(x) = (x + |x|)/2 is the ramp function, V is a differential operator: 

V := Z x dy 2 d Ys + Z 1 Z 2 dy 3 d G + cyclic + Z x Z 2 Z 3 , 
and ££i is the lattice: 

&i ■= { (n, 75, r 3 ) E N 3 I r i+ i + Ti_i < A*i } ■ 
(P)M - (P)M - Massive 



(1.13) 



(1.14) 

(1.15) 
(1.16) 



For arbitrary ^3 : 



Ml M2 



C = E ^ ^ ^ 1CT2 (^3, #1, #2, # 3 ) 

<Tl=0 (T2=0 



(1.17) 



For Hi + — Hi-i S 2Z, one also has 



c= z " z ? z " Y l 

(ri,r2,r 3 )eifin^2n^3 

xe~i v K T1T2T3 (Y,G) 
(P)M-(P)M-(P)M 
For arbitrary ^'s : 



R{r 1+ n -^-^ ) a^+ gaziazgi ) r(t 3 + ^-^-^ ) 

(1.18) 



G=G(Y,Z) 



Ml M2 M3 



(Tl =0 CT2 =0 (T3 =0 



(1.19) 



For /ij + Hi+i — Hi-i 6 2Z, one also has 



c= z i T1 z ? z ? y i 

(ri,r2,r 3 )eJ^inif 2 nj% 

xe"i p r iT2T3 (7,G) 



fl(n+ ^-^-^ ) r(t 2 + ^-m3-m^ yfl ^ + >m ) 

Yc 



G=G(Y,Z) 



Y, 



(1.20) 
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Organization of the paper 

Section 2 contains some general discussions about polynomial solutions to PDEs. In Section 
3, we provide the most general solutions to one PDE encoding the interactions of one (P)M 
field. The solutions to the system of PDEs corresponding to general couplings that involve 
more than one (P)M field are derived in Section 4. Our results are discussed in Section 5. 
Finally, Appendices A, B, C and D contain some mathematical details on the derivations 
of the results presented in this paper. 

2 General structure of the solutions 

Before solving the consistency equation (1.8), let us first discuss in some detail the general 
idea underlying our way of organizing the polynomial solutions. The latter are in fact the 
only relevant solutions for the analysis of HS cubic interactions. 

2.1 Cubic interactions as a vector space 

Restricting the attention to polynomials in Y and Z , the function C(Y, Z) can be expanded 
as 

C(Y,Z)= C^ a Z?Z?Z?Y^Y^Yp, (2.1) 

<Ji,Ti>0 

where C£*p^)? are arbitrary coefficients. Since the operators d (1.9) preserve the spin 
degrees: 

Si = Oi + Ti+l + Ti_l , (2.2) 

any solution to eq. (1.8) can be decomposed into solutions with fixed spins Sj . Moreover, 
being a linear PDE, for given Sj its solution space is a finite dimensional vector space, 
whose elements can be labeled as: 

N 

C s lS2S3 ^ Z ) = Y J K n P° 1S2S3 (Y, Z) [K n G R] . (2.3) 

n=l 

Here, N is the dimension of the vector space, namely the number of s\ — S2 — S3 independent 
couplings, whereas 

{ Pl 1S2S3 (Y, Z) , P 2 SlS2S3 (y, Z), ... , P S N 1S2S3 (Y, Z) } , (2.4) 

is a corresponding basis. Therefore, solving eq. (1.8) is tantamount to finding a set of basis 
vectors P^ lS2S3 (Y, Z) . Let us notice that, since the spins do not enter explicitly the PDE, 
its solutions are spin independent. Hence, one can first determine the basis {P n (Y, Z)} 
without specifying its spin dependence, and then restrict the attention to fixed spins. 

2.2 Choice of basis 

The basis {P n (Y,Z) } is of course not unique, and some choices can be more convenient 
than others. In order to clarify this point, let us first expand the basis solution P n (Y, Z) 
in powers of 8/L : 

P n [i;Y,Z)=Y,(i) k P^(Y,Z)- (2-5) 

k>0 



- 5 - 



In physical term, this is an expansion in the number of (ambient-space) derivatives, so that 
the leading term P n 0) can be identified with the highest-derivative piece of the coupling 
P n (Y, Z) . A convenient choice is a basis whose elements P n (Y, Z) have all different leading 
terms: 

P^(Y,Z)^P^(Y,Z) Wn^m. (2.6) 

Let us notice that this requirement does not fix the basis uniquely. Indeed, denoting by A n 
the maximum number of derivatives in P n (Y,Z), one can always add a lower-derivative 
piece P m (Y, Z) as 

/ » \ |(A n — A m ) 

Q n (Y,Z) = P n (Y,Z)+^y P m (Y,Z) [A n >A m ], (2.7) 

without spoiling the condition (2.6). Our strategy in order to construct a basis of this type 
consists in two steps: 

• identify all possible leading terms Pn (Y,Z); 

• for each P n 0) (Y, Z) , find (if it exists) a corresponding full solution P n (Y, Z) . 

As we have just explained, the solution P n (Y, Z) is not unique but can be always replaced 
by some other solutions Q n (Y, Z) differing by lower-derivative solutions. In the next Section 
we exploit this freedom in the choice of basis in order to simplify our analysis. In particular, 
we first construct a basis: 

B P = {P n (Y,Z)}, (2.8) 

proving that it spans the entire solution space of eq. (1.8), and then we introduce another 
set of solutions: 

B Q = {Q n (Y,Z)}, (2.9) 

whose elements satisfy Qn\Y, Z) = Pn\Y, Z) for all n 7 (so, the index n can be understood 
as a label for different leading terms). The latter conditions ensure that B Q be also a basis 
of the solution space, i.e., Span(B Q ) = Span(l?p) . The reason for this change of basis is 
that, while B P is convenient to prove the completeness of the solution space, B Q turns out 
to be more suitable in constructing solutions to more than one equation. 

3 Solutions to one equation 
3.1 Massless equation 

The aim of this Section is to find the general solution to eq. (1.8) . Since the latter consists 
in a product of commuting operators L\ , it is convenient to first analyze the kernel of a 
single operator. In particular, we start analyzing the massless case (fi\ = 0), where eq. (1.8) 
reduces to 

C 1 (fi 1 )C(X,Z) = 0, (3.1) 

with the operator C\ given by 

A (Ai) = y 2 d Zs - y 3 o Z2 + 1 (y 2 d Y2 -Y 3 dy 3 + f ) d Yl . (3.2) 

7 More precisely, we shall provide linear combinations satisfying this equality. 
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3.1.1 General solution 

Eq. (3.1) can be solved as a power series in Y\ : 



C n (Y, Z) = Y J C%(Y 2l Y 3 ,Z){-±d Yl ) k Y^ , 



(3.3) 



k=0 



where C cri (1", Z) denote the solutions whose highest power of Y\ is o\ . Plugging the 
expansion (3.3) into eq. (3.1), one ends up with a differential recurrence relation for C< 



(fc) . 



(Y 2 d z , - Y 3 d Z2 ) Cg>(Y 2 , Y 3 , Z) = (Y 2 dy a -Y 3 dy 3 + % /2l) C^\Y 2 , Y 3 , Z) , (3.4) 

where C^ n 1) = . The latter can be solved iteratively starting from the leading term: 

Cg(Y 2 ,Y 3 ,Z) = K{Y 2 ,Y 3 ,Z l ,G 1 {Y,Z)) 

= ^ W ^ ly / 2 ^ 3 [Gi(y,2)f, (3.5) 



where Gi(Y, Z) :=Y 2 Z 2 + Y 3 Z 3 . After summing C ai over o\ , the full solution reads (see 
Appendix A for details): 



C(Y,Z) 



E 



K 



(3.6) 



01 ,<J2 ,(T3 ,71 ,1) 



where the polynomial functions Pa\a 2 a z v are defined as 



Pala 2 a 3 v(^l'i Y, Z) — 



E 



Z T! jmti y^a j^oa [ Gl ^r ^ 



(3.7) 



and [o] n := a(a — 1) • • • (a — n + 1) is the descending Pochhammer symbol. Let us empha- 
size that, due to the presence of the terms [a 2 + cr 3 ] p+q in the denominator, these functions 
may be ill-defined. More precisely, potentially diverging terms appear for those values of 
p and q such that 

(?2 + o 3 + 1 < p + q < o\ ■ (3.8) 
Hence, there are only two cases in which the solutions are well-defined: 

• When a 2 + a 3 > u\ . In this case, the problematic terms are simply absent and one 
gets the following set of solutions: 



BpiGh) = { P£ niraV (fii) | <r 2 +a 3 >a 1 } 
• When the condition (3.8) holds and the residue: 



(3.9) 



h + fL>3-f] ? = (f + <r 2 )---(f -a 3 )x 

x (-!)<? (f + q - a 3 + 1) • • • (f - p + a 2 + 1) , (3.10) 
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of the corresponding diverging piece vanishes. This happens when 

(f + <7 2 )---(f -<x 3 )=0, (3.11) 
so that the following set of solutions: 

Bptifii) = { P^ nnv {fli) | ^3 > f , f S Z> ; a 2 > |f | , f e Z <0 } , (3.12) 
are well-defined. 

Notice that Bp 2 is non-empty only when \x\ is an even integer. Moreover, it has a non- 
vanishing intersection with Bp\ , so much so that, when Jx\ = , Bp\ becomes a subset of 
Bp2 ■ Finally, the union of the two sets : 

B P (fl 1 )=Bp 1 (j2 1 )UBp 2 (j2 1 ), (3.13) 

forms a (redundant) basis 8 of the solution space of eq. (3.1). Comparing eq. (3.13) to 
eq. (2.8), one can see that the index n in P n (Y, Z) corresponds to the collective index 
Wi, &2, 03, ti, v} labeling different leading terms. 

3.1.2 Change of basis 

At this stage, we have explicitly constructed the basis (3.13) of the solution space of 
eq. (3.1). However, although the form (3.7) in which the basis solutions are written makes 
the completeness of the solution space manifest, it is not suitable for the analysis of more 
than one equation. Therefore, in the following we construct other two sets of solutions 
Bq and Bq , and, analyzing their leading terms, we prove that their union spans the same 
space as Bp . 

Q solutions Let us notice that the operators: 9 

Hi := d Xi+1 ■ d x ^ d Ui+1 ■ du^ - • d Xi+1 d Ui+1 ■ d x ^ , (3.14) 

commute with the gradient operators: 

[H i ,U r d x .} = 0, (3.15) 

without relying on the on-shell conditions. As a consequence, one can easily construct 
couplings of the form: 

K(Y 2 ,Y 3 ,Z 1 ,H 2 ,H 3 ), (3.16) 

that are invariant under the gauge transformations associated to the first field. These 
couplings are not written as functions of Y and Z , but they can be always brought to that 
form by performing the integrations by parts of all the total-derivative terms present in H 2 
and H 3 . Hence, one can consider the following set: 

Bq{H) = {Ql\„ zh2h M)} > ( 3 - 17 ) 

8 With a slightly abuse of notation we use the term 'basis' to denote a set of solutions which spans the 
entire solution space, regardless of the fact that they be linearly independent or not. 

9 These operators can be obtained as deformations of their fiat-space counterparts (see Appendix D). 
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whose elements are given by 

Q T a M M^ Z ) ~ Z " Y P *2 2 H 3 k3 , (3.18) 

where ~ means equivalence modulo integrations by parts. Although finding the exact form 
re< l u ^ res an involved analysis (see Appendix C for details), for our purpose it 
is sufficient to identify the corresponding leading terms: 

Q^tohsiMi; Y, Z) = Z? Y^ +h * Yp +f * Y°* +h > + 0[l) ■ ( 3 - 19 ) 

At this point, it is straightforward to check that these leading terms exactly reproduce the 
ones of P£l a2(T3V in Bpx (3.9) for v = 0. The dependence on G\ can be recovered by taking 
proper combinations of them. For instance, the 3 — 2 — 3 coupling starting with the leading 
term Y\ Y 3 G\ is given by 



^ (*2 h 2 ) 2 - /ii (y 2 H 2 ) (Y 3 H 3 ) + (y 3 H 3 f 



Ho 



- (!) 2 ^ f ^ [ y i y 3 + | z 2 Gi y 2 ^ 2 + ^ y 3 ^)] • (3.20) 

As one can see from this example, depending on the values of \1\ , the relation between the 
elements of the two basis Bq(Ji\) and Bp\{fix) can be singular. However, the vanishing 
coefficient (i.e., (pi + 2) \1\ (fix — 2) in (3.20)) is overall so that one can always normalize the 
leading term. See Appendix C for the general case. Since all the leading terms of Bpi(jli) 
can be reproduced by Q^ 2a3 h 2 h 3 (P'^) > an< ^ ^pi(Pi) covers the whole solution space when 
(1.10) is not satisfied, one can conclude that 10 

Span(5^(pi)) = Span(£pi(/2i)) . (3.21) 

Q solutions Let us notice that, when the constant in eq. (3.1) is an integer number, 
it can be 
function 



it can be removed factoring terms of the form Y 2 Ml ^ 2 or Y^ 1 ^ 2 . Hence, using the ramp 



R(x) := (\x\+x)/2, (3.22) 
one can construct the following solutions: 

C(Y,Z) = Y 2 R{ ^ Y 3 R ^ e'T^ K(Y 1 ,Y 2 ,Y 3 ,Z 1 ,G) [if f eZ], (3.23) 

G=G(Y,Z) 

with 

G(Y, Z) := Yi Z x + Y 2 Z 2 + Y 3 Z 3 . (3.24) 

Here e~L V K are the general solutions to the massless equation £i(0) C = , where the 
operator T> is defined as 

V := Zx dy 2 dy z + Z x Z 2 dy 3 d G + cyclic + Z x Z 2 Z 3 8% . (3.25) 



10 By continuity of the solution spaces Bpi(fj,i) and Bq([Ai) in fii , this statement holds also when the 
condition (1.10) is satisfied. 
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The solutions (3.23) can be decomposed in terms of the following functions: 



G=G(Y,Z) 



= Z? Y* 1 Y° 2+R ( ^ Yp +R ^ ) [Gi (Y, Z) + Yi Z x f + , (3.26) 

whose leading terms coincide with those of P£l a2(T3V in -Bp2 (3.12). As a consequence, the 
set: 

B Q (fii) = {Ql\ a2CTs M}, (3-27) 
together with £?q(Mi) span the entire solution space: 

Span^Gux) U Bqfjix)) = Span(-B P i(/Ii) U £ P2 (Mi)) • (3.28) 
3.2 Partially-massless equation 

Let us now move to the case in which one PM field (/ii € N) is involved in the interactions. 
Then, the corresponding cubic vertices have to satisfy the PDE: 

A (Ml " Ml) A (Ml - Ml + 2) • • • A(Mi + Mi) C(y, Z) = . (3.29) 
3.2.1 General solutions 

The general solutions to the above equation can be decomposed in terms of the functions 
CVs satisfying 

C\{fii-Hi + 2k)C k = 0, (3.30) 

with k = 0, 1, . . . , n\. The CVs are given by (3.6) where fii is replaced by Mi — Mi + 2 k , 
hence the following set of functions: 

{ Pala 2 a 3 vifil ~ Ml) } U { P^ 2a3V (fil ~ Ml + 2) } U • • • U { P^.^l + W ) } , (3.31) 

forms a complete basis for the solution space of eq. (3.29). However, this decomposition 
is possible only when o"2 and (73 are regarded as real (non-integer) numbers. Indeed, as 
in the massless case, the set (3.31) contains solutions that are ill-defined. Therefore, out 
of them, one needs to select only the ones that are well-behaved as o<i and 03 approach 
some integer numbers. In the massless case, for any given leading term there is only 
one function in the set, so that the selection simply amounts in examining whether each 
function Pala 2 <j 3 v(P'i) is well-defined or not. On the other hand, in the PM case one has 
H\ + 1 solutions -PJJ 0-20-3^(^1 — Mi + 2 k) with the same leading term, so one has to analyze 
all possible linear combinations of them. 

PI solutions Let us consider the first divergent term in the series expansion of Pa[a 2 viv( x ) 
(3.7). It arises for p + q = a% + 0-3 + 1 and its residue is (| + 02) ■•■(§ — 03) (— l) 9 (see 
eq. (3.10)). Therefore, it can be cancelled by taking the following linear combination: 

(f + 1 + a 2 ) P T a \ a2a3V {x) - (f - a 3 ) P^ 2(T3V {x + 2) . (3.32) 
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Although one can get rid of the first divergent term, there is no way to eliminate the 
remaining divergent pieces that are still present for p + q > 02 + 03 + 2 . Hence, the com- 
bination (3.32) is well-defined only when the order of the polynomial is low enough in order 
to prevent those divergent terms from showing up, that is, when 02 + 03 + 1 > a\. Let us 
notice that this condition on the leading term is very similar to the one obtained in the 
massless case (3.9), but for the shift by one on the left-hand side of the inequality. In 
general, taking the following linear combination: 

pN r i ( r \ ■ = 

= E [l + n + ^ln- k (f " K\„ 2 „A* + 2 *0 > (3-33) 

one can cancel the first n divergent terms. Once again, as these functions contain divergent 
terms for p + q > a 2 + 03 + n + 1, only the ones belonging to the set: 

B [ ;\{x) = { PW£ a3V (x) I a 2 + a 3 + n > a x > n } , (3.34) 

are well-defined. Here, the lower bound on 01 has been introduced since the solutions 
Paia2(T3v with 01 < n — not having a sufficient number of terms — are simply linear combi- 
nations of PaiJ^a 3 v , PaiJaazv 1 ■ ■ ■ 1 P^la^llv ■ To sum up, starting from the formal basis 
(3.31), one can extract the following set of well-defined solutions: 

:= B [ p[(fll ~ Ml) U B$[(fii - Ml) U ■ ■ ■ U - Ml) • (3.35) 

P2 solutions On the other hand, precisely as in the massless case (3.12), the divergent 
terms are harmless when the corresponding residues vanish. Since everything is well-defined 
in this case, the corresponding solution space is spanned by the union of the massless sets 
(3.12): 

B P2 {fii - Mi) U B P2 (fl 1 - mi + 2) U • • • U £ P2 (mi + Mi) • (3-36) 

However, for the succeeding analysis it proves convenient to change basis. For this purpose, 
let us first introduce the new basis functions: 

fc=o ^ ' 

which are linear combinations (with binomial coefficients) of \n\ + 1 functions with con- 
secutive arguments starting from x . The functions (3.37) are well-defined for n > if 
03 > x/2 + n and for n < if 02 > x/2 — n . As we will see in the next Section, these lin- 
ear combinations allow the compensation of the first \n\ terms in the Yi-expansion, making 
the link to the PM counterpart of the Q basis function more straightforward. In terms of 
these functions, one can define the following sets: 



Bp2(mi,Mi) 



Paiwsv (-R(Mi - Ml)) I °3 > n > [n>0\ 

, (3-38) 

<Paia 2 (T3V (-i?(-Mi -Ml)) I °2 > \n\ j [n < OJ 
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such that their union: 



Bp2(m,fii) :=B p2 2 U B p2 



J (/xi,/xi) U ••• U B 



P2 



(3.39) 



spans the same space as (3.36). To be concrete, let us consider the //i = 2 case, where one 
can distinguish between three different subcases: u\ > 2 , p,\ < — 2 and /Zi = . If /Zi > 2 , 
then one gets 



r Mi-2 i 

R l 2 1 

n P2 
B P2 



{ P^ T1 ffi, 

L}n 



2) = K 



<7 3 > 



(Ai - 2) 

{^CTi(T20"3»j(Al — 2) = -PcrJeraiTalXMl — 2) — P a \ a2a ^ v (fil , 



> 



Ml \ 
2 J » 



(3.40) 



(mi - 2) - 2P£ f/i!) + 



0"10"20"3t 



CI (T2 0"3U 



(Ml + 2) 



while, if /^i < —2 the corresponding sets are given by 



r Mi-2 i 

M 2 



'P2 



B l 2 



pri 

<T1CT2(T3^ 



(Al + 2)-2i^ oaa8V (Ai) + i^ oa , T , 



u(Ml 



£72 > 



»i-2|l 
2 1/ ' 



F2 



{Pt<S£v(fil + 2) = P^vifil + 2) - P? 1<72 a 3V (fil) | ^2 > if |} , 



B P2 2 



{pW^J/il + 2) = P? ia2 „ 3V (fii + 2) 
Finally, in the last case R\ = 0, one ends up with 



} 



(3.41) 



4f = {4^fe(o) | *2 > i} , 41 = {^ fc 3 ,(o)} > *K = 



^3 > l} • 

(3.42) 

As one can see from the above example, the set (3.38) is chosen such that its elements are 
linear combinations involving always the P^ a2(T3V (x) with x = m.m{\u\ — . . . , |/2i + 

All in all, the solution space of the partially-massless equation (3.29) is spanned by 
the union of the two sets (3.35) and (3.39). 

3.2.2 Change of basis 

After having generalized the sets Bp\ and Bp 2 to the sets (3.35) and (3.39), we now aim 
at finding the analogue of the convenient basis Bq (3.17) and Bq (3.27) for the partially- 
massless case. 

Q solutions The first observation is that the identity: 
[Y 1 n K(Y 2 ,Y 3 ,Z 1 ,H 2 ,H 3 ), (C/i-a Xi r +1 ] =0 



n = 0,l,..., W ], (3.43) 



holds because the operator K(Y2, Y3, Z%, H%, H3) commutes with U\ ■ dx 1 (without on-shell 
conditions). Thus, let us consider the functions: 



(3.44) 
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that are of the same form as (3.18) but for an additional factor of Y™ . In the previous 
sections, we have shown that, for n = 0, the sets Bq (3.17) reproduce all the leading terms 
with <72 + 03 > o\ in Bp\ . For n > 1, it is still manifest that the following set: 

gives leading terms with a 2 + 03 + n > o\ > n and v = , reproducing the same bound 
(3.34) on u\ . However, for a given leading term, depending on the values of the <7j's, there 
can be more than one solutions — more than one values of n can satisfy the inequality — in 
both (PI and Q solutions) sides. In order to handle this subtlety as well as the v > 1 
case, one can consider once again a change of basis in such a way to remove the degeneracy 
between the leading terms. Concerning the PI solutions, this can be done by taking 
particular linear combinations of the Pj-^ojrsw's, with n = fix, fix — 1, . . . , y«i — m , that give 
rise to the cancelation of the first m leading terms. A similar analysis can be also done for 
the Q solutions after integrating by parts the total-derivative terms present in H2 and H 3 
(see Appendix C for more details). In the end, one can see that all the leading terms of 
the PI solutions in the set (3.34) can be reproduced by the Q solutions in the sets (3.45). 
Hence, exploiting the completeness of the PI solutions when the condition (1.10) is not 
satisfied, one can conclude that 

Span [.Bq (/Zj. , //1 )] = Span[Bpi(/ii,/ii)] , (3.46) 

where 

Mi 

B (/2 1)Ml ):= U^'fe.w)- (3-47) 

n=0 

Q solutions Here, we search for a generalization of the solutions (3.23) to the partially- 
massless case. 11 When ^1 = 1, due to the identity: 

- 1) + 1) z 2 = [Z2 A (Pi + 1) - 2 Y 3 ] £i(/ii - 1) , 
d {Jii -l) A (Ml + 1) Z 3 = [Z 3 d (fix - 1) + 2 Y 2 ] + (3.49) 

one can easily check that the functions: 



C(Y,Z) = Z 2 Y 2 K 2 J Y 3 V 2 > e-T v ' Ki(Y 1 ,Y 2 ,Y 3 ,Z 1 ,G) 



+ Z 3 y 2 v y y 3 v ' e-T u K 2 (Y 1 ,Y 2 ,Y 3 ,Z 1 ,G) , (3.50) 

G=G(Y,Z) 

solve the PDE £i(/2i — 1) Ci(fii + 1) C(Y, Z) = . In general, one can consider the set: 

P^ T3l (/2 1 ,Mi) = {Ql^}> (3-51) 



11 Because of (3.28) and (3.36), one is naturally led to consider the basis: 

B Q (p,i - mi) U Bq(mi - Mi +2) U ••• U BQipi + m) . (3.48) 
However, the latter is not suitable to study the solution space of more than one PDE. 
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whose elements are given by 



r>l T 2> T 3] T i (V 7\ — 7 r i 7 r 2 7 T 3 V V 2 2 
Vaicra^^' Z ' ~~ Z l Z 2 X 3 r 2 



xe-L V Y 1 ai Y 2 a2 Y 3 CT3 G v 



(3.52) 

G=G(y,Z) 

with T2 + T3 < Ml (see Appendix B for details). The proof that the union of these sets 
together the Q solutions span the entire solution space relies once again on the leading 
term analysis. Using the identity (A. 9), one can show that 

P^ 3n (x) = (l) H [(Z 2 d Y3 -Z 3 dY 2 ) H Z?Y?Y**Y°* [G 1 (Y,Z)Y' + 0{1) 

= (!) H e (-ip (';') N.a x 

T2+T 3 = \n\ V 3/ 

x Z? Z? Zp Y* 1 Yp~ Ta Y^ T2 [Gi (Y, Z)] v + O ( (|) H+1 ) . (3.53) 

In this form, it is straightforward to check that each leading term above can be reproduced 
by the functions (3.52). As an example, let us assume Mi > Ml . In this case, the function 

(3.53) in B p2 2 (mi,Mi) has to satisfy the condition 03 > (mi — Ml)/2 + n. As a con- 
sequence, the minimum power of I3 is (Mi — Mi)/2 + T3 , coinciding with the one in (3.52). 
All in all, we have 

Span -Bpi(Mi,Ml) U B P2 (p.i, Ml)] = Span ^Bg (pi, Mi) U Bq(Mi,Mi) , (3.54) 

where 

B g (/ii,/ii):= |J Sg 8 ' 781 ^!,^). (3-55) 

T"2+T3</il 

Let us mention that, though highly redundant, the basis (3.55) proves very convenient in 
the study of the solution space associated to more than one PDE. 

Before moving to the next Section, let us summarize the general solutions to eq. (1.8) 
in the Q and Q basis. 



Solutions to one equation 

For arbitrary M2 — M3 : 



C = E ^ ^ ( y2 > y 3> Z X ,H 2 ,H^) 
0-1=0 

For Mi + H2 — M3 6 2 Z , one also has 



(3.56) 



C 



\ 7T1 7^2 V V 2 / V" V / 

2 , ^1 Z 2 ^3 r 2 r 3 



(Tl,T2,T3)e-S?l 



G=G(y,Z) 



where 



x e -i v K T1T2T3 (Y,G) 
££i := { (7-1,75, r 3 ) E N 3 I r i+ i + n-i < Mi } • 



(3.57) 
(3.58) 
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4 Intersection of the solution spaces 

Let us now consider general interactions involving more than one (P)M field. In these 
cases, depending on their number, one has to consider the intersections of the solution 
spaces of the corresponding PDEs. In the following, we carry out this analysis for the Q 
and Q solutions separately. 

4.1 Q solutions 

When two (P)M fields (say i = 1,2) are present, one has to take the intersection between 
the coupling (3.56) and its cyclic permutation: 

Ml M2 

Y* 1 (Y 2 , Y 3 , Z X ,H 2 ,H 3 ) , Yl Y 2 2 k ? (Y 3 ,Y 1 ,Z 2l H 3 , H x ) , (4.1) 

ui=0 U2=0 

that is 

Ml M2 

C = E ^ ^ ^ CT1<T2 (Y 3 , H X ,H 2 , H 3 ) . (4.2) 

(Tl=0 O"2=0 

When all three interacting fields are (P)M, one has to intersect the couplings (4.2) with 

M3 

]T Y^K^(Y 1 ,Y 2 ,Z 3 ,H 1 ,H 2 ), (4.3) 

o-3=0 

leading to 

Ml M2 M3 

C = ^ J2 ^ Y 2 2 Y z 3 K aia2a:i (H l7 H 2 ,H 3 ) . (4.4) 

CTl=0 CT2=0 CT3=0 

4.2 Q solutions 

For the analysis of the intersections among the Q solutions (3.57), it is convenient to start 
with the leading terms C (0) of the solutions to the entire system of PDEs. Plugging the 
expansion (2.5) into eq. (1.8), one discovers that C (0) satisfies the following relatively simple 
equations: 

{Yi+i d Zi _, - y,_x d Zi+1 f* +1 C<°>(Y, Z) = 0, (4.5) 
which can be solved in terms of an arbitrary function K T1T2T3 (Y, G) as 

C {0) (Y, Z)= Z i T1 Z 2 2 z z K T1T2T3 (y, G{Y, Z)) . (4.6) 

(ri,T 2 ,T3)ei? 

Here, the lattice Jz? is given by the intersection of the lattices Jzfj's (3.58) associated with the 
i-th PM fields. In the following, we derive the general solutions to the system of eqs. (1.8) 
starting from the leading terms C (0) of eq. (4.6) . For that, we distinguish between the 
cases where two or three (P)M fields are involved. 

When two (P)M fields (say i = 1, 2) are present, one has to solve the system: 

A(M3 - M2 - Mi) A(M3 -M2 - Ml + 2) • • • A(M3 - M2 + Mi) C(Y,Z) = 0, (4.7) 
£ 2 (mi - M3 - M2) £ 2 (mi - M3 - M2 + 2) • • • £ 2 (mi - M3 + M2) C(Y, Z) = , (4.8) 
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where the leading terms of the corresponding couplings are given by (4.6) with 

if = jgfi n i? 2 • (4.9) 

At this point, it is convenient to split the sum over (ti,T2,T3) into the two regions: 

^in^ 2 ni?j, ^ni^n^, (4.10) 

where J2?j C denotes the complement of the set Jzfj . When (n, T2, T3) € J2?i n ^2 H J2?3 , out 
of the three combinations Tj + ^i-v-i+i-Hi-i ^ a ^ mos ^ one can De positive. Hence, without 
loss of generality, one can assume 



Then, depending on the sign of n + 2 ~ 5 there are two subcases. When the latter is 
non-positive, a comparison with eq. (3.57) and its cyclic permutation directly shows that 
the coupling: 



C(Y, Z) = Z{ 1 Z? Zp e-T v K T1T2TS (Y, G) 
is a solution to both PDEs (4.7,4.8). On the other hand, when 



G=G(Y,Z) 



one can consider the following function: 



C(Y, Z) = Z? Z 2 T2 Z^ Y 1 



e -L u K TlT2T *(Y,G) 



G=G(Y,Z) 



(4.12) 



(4.13) 



(4.14) 



that is a solution to eq. (4.8). Although it is not manifest, the latter solves also eq. (4.7). 
One way of proving it consists in pushing the exponential function to the left, ending up 
with: 



C(Y,Z)= £ Z?Z?Zpe~LVK^p(Y,G) 

p,g,r>0 

where fi = T2 + p + r , T3 = T3 + q + r and 



G=G(Y,Z) 



(4.15) 



p 9 r a y 9 a/5 G r K nr2T3 (y,G). 



(4.16) 



Comparing with (3.57), one can see that the latter is a solution to the first equation 
provided the conditions: 



n + Tl < A*3 , t 2 + r 3 < /ii , r 3 + n < /i 2 , 
f 2 + fQrfarfQ < , r 3 + < ? 



(4.17) 
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are satisfied. These conditions hold for any (t%, T2, T3) G Jz?i fl J*?2 H satisfying conditions 

(4.11) and (4.13), and therefore the function (4.14) solves eq. (4.7). Finally, the solutions 

(4.12) and (4.14) can be written at once in a cyclic form as 



C 



Z? Z? Z? Y, 




2 




2 




M3-M1-M2 
2 



) 



X 



(Ti,T 2 ,T 3 )eifirL%n.% 



x e i 



K TlT2T3 (y,G) 



(4.18) 



G=G(Y,Z) 



Let us now consider the case (ri , T2, T3) € Jz?i fl ^2 H Jz? 3 c , where ^3 < ti + T2 < S3 . If one 
requires that the third field be a unitary massive field, then this type of interactions are 
ruled out in dS (see Figure l). 12 On the other hand, in AdS these would correspond to the 
interactions between two massless and one massive fields with /1/3 < , but we do not find 
any solution of this type. 

Last of all, when all three fields are (P)M, the leading terms of the corresponding 
couplings are given by (4.6) with 



This situation is nothing but a subcase of the two (P)M interactions that we have considered 
before, therefore the solution is given by (4.18). 

At this point, we have completed the analysis of the (P)M cubic interactions, and the 
results are summarized in Section 1.2. 

5 Discussions 

To conclude, we discuss the implications of the condition (1.10) in more details, pointing 
out the key differences with respect to the flat-space case. Moreover, in order to make 
contact with our previous work [35], we provide an example. 

Non-Abelian interactions Let us recall that the cubic interactions which exist regard- 
less of the condition (1.10) are given by arbitrary functions of H (Q solutions). Since the 
latter are trivially gauge invariant — their gauge invariance does not rely on the on-shell 
conditions — they do not lead to any deformation of the gauge transformations. To iterate, 
they are all Abelian and of the Born-Infeld type, expressible in terms of linear curvatures. 
When the condition (1.10) is satisfied, besides these //-couplings, supplemental G-couplings 
(Q solutions) appear. However, only a part of them is independent from the //-couplings 
and may lead to non- Abelian deformations of the gauge symmetries (see Figure 2 for a 
schematic picture). Let us stress that there are many PM interactions which do not sat- 
isfy the condition (1.10), including notably the PM spin-2 self-interactions [61, 63]. Those 
interactions cannot give rise to any non- Abelian deformations of the gauge symmetries. 

12 The third field may be PM, but then one has to impose gauge invariance also under its gauge symmetry. 
This case is considered later. 



if = i # 1 n^ni , 3. 



(4.19) 
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Figure 2. H- and G-couplings. Non-abelian couplings are a subset of the G-couplings that can 
not be written as -ff-couplings. 

Interactions of one massless and two massive fields Apart from PM interactions, 
which do not have any flat-space counterpart, a novel property of (A)dS is present in the 
cubic interactions of one massless (fj,% = 0) and two massive (/j, 3 > ^2) fields. According 
to the condition (1.10), supplemental G solutions appear when ^ - /J2 £ 2Z> . Among 
them, the non- .ff-couplings are the ones satisfying conditions (3.8) and (3.12): 

OX > (T 2 + (7 3 + l, C7 3 >i^, (5.1) 

where the <7j's are related to the spins Sj's as 

Sl = <T\ + V , S 2 = a 2 + V + Tl , S3 = CT 3 + V + T\ , (5.2) 

with f and t\ being the powers of Gi and Z\ respectively (see (3.7)). This type of interac- 
tions is related to non-trivial Noether currents. In particular, the electromagnetic (1 — s—s) 
and the gravitational (2 — s — s) minimal couplings correspond to \i2 = A*3 ; as in flat space. 
On the other hand, there are many interactions that satisfy the above conditions for pos- 
itive even integers /is — fj,2 ■ For instance, when /X3 — \x<i = 2 , one can have 2 — s — (s + 1) 
interactions involving non-trivial Noether currents made of fields with different masses. 
This is a novelty with respect to flat space, where the interactions leading to non-trivial 
conserved currents are only available when the two masses are equal: ^2 = ^3- This obser- 
vation suggests that, in (A)dS, (HS) multiplets may involve not only fields with different 
spins, but also with different masses. Let us also mention that the minimal-like couplings 
si— s — s require /X3 — ^2 < s\ — 1 , so that in this case interactions involving fields with 
different masses are available only for si > 3 . 

Example: 4—4—2 interactions For concreteness, starting from the general solutions 
provided in this paper, we show how to recover the example of 4—4—2 interactions obtained 
in [35] by means of a numerical algorithm. The latter are the cubic interactions between 
two spin-4 fields at their first PM points = \i% = 1) and a massless spin-2 field (/is = 0). 
In this case one can find two Q solutions (1.19): 

H x H 2 Hi, H 1 H 2 HiY 1 Y 2 , (5.3) 

and six Q solutions (1.20): 

Z 3 e-i v Y*Y 2 3 Y 3 2 , 
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■i: v Y 1 d Y 2 6 Y 3 G, 



l u Y{ Y 2 Z G , 



Z 3 e-t v Yi 2 Y 3 G, Z 3 e-i u Y x Y 2 G 



v - 



(5.4) 



However, two of the above Q couplings — the first and the fourth — can be expressed in 
terms of the Q couplings (5.3) and the remaining Q couplings. Hence, one is left with six 
independent solutions (two -ff-couplings and four (non-H) G-couplings) in complete agree- 
ment with the result of [35]. More explicitly, for instance, the lowest-derivative interaction 
Z 3 e-T V Y 1 Y 2 G 2 gives 

Z 3 Y X Y 2 G 2 -{ Z\ [G 2 + 2 {Yx Z x + Y 2 Z 2 ) G + 2 Y x Y 2 Z x Z 2 ] + 4 {{f Z x Z 2 Z 3 3 , (5.5) 

reproducing the vertex in eq. (3.40) of [35]. 

Acknowledgments 

We are grateful to A. Sagnotti and E. Skvortsov for discussions. The present research 
was supported in part by Scuola Normale Superiore, by INFN, and by the MIUR-PRIN 
contract 2009-KHZKRX. The work of MT was also supported by the Max Planck Institute 
for Gravitational Physics (Albert Einstein Institute). 

A Recurrence relation 

In this Appendix, we provide some details on the derivation of the solutions to the massless 
eq. (3.1) and to the PM eq. (3.29). 

A.l Massless equation 

As explained in Section 3, plugging the expansion: 

oo 

C ai (Y,Z) = J2c^(Y 2 ,Y 3 ,Z) {-idy^Y^, (A.l) 

k=0 

into the equation: 

Y 2 d z , - Y 3 d Z2 + | (y 2 d Y2 -Y 3 dy 3 + f ) %] C(Y, Z) = 0, (A.2) 
one ends up with a differential recurrence relation for : 

(Y 2 d z , - Y 3 d Z2 ) C£(Y 2 , Y 3 , Z) = (Y 2 dy 2 - Y 3 d Y:i + \ fix) C^(Y 2 ,Y 3 , Z) , (A.3) 
where C^ 1 ' = . The latter can be solved iteratively starting from k = : 

C^(Y 2 ,Y 3 ,Z) = K(Y 2 ,Y 3 ,Zx,Gx(Y,Z)) 

= E K Zwsv Zi 1 Y2 2 Y-P [Gi(Y Z)Y , (A.4) 

where G\(Y, Z) := Y 2 Z 2 + Y 3 Z 3 . For this purpose, one can consider the following ansatz: 

C ai (Y, Z) = J2 Klwsv Pala 2 a 3 v(fii;y, Z) , (A.5) 
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with 



P^ 3V (jii;Y 2 ,Y 3 , Z) = J2 cti (Z 3 dY 2 ) k ~ e (Z 2 dy z Y Y°* Y.p [Gi(Y, Z)f , (A.6) 



and turn eq. (A. 3) into a recurrence relation for c k £ ■ 

(k - I) C k; £ + {I + 1) C k) i+1 = C k -1,£ , 

(k - £) (<t 3 - £) c k , -(£+l)(a 2 -k + £ + l) c fc 



Mi „ 

— Ck-1,£ ■ 



The latter can be straightforwardly solved as 

1 [ CT 2+ 2 J/,-, 



(A.7) 



(A.8) 



(k-t)H\ [*2 + <T3} k 

Plugging the coefficients (A.8) into eq. (A.6) and summing over k, one gets eq. (3.7). 

To conclude this part of the Appendix, let us provide an identity involving the function 

pri 



Pj-i (*-n + 2fc) 



|) n (Z 2 dy 3 -Z 3 dy 2 y 



X 



E 



iz 3 d Yl d- 



Y-2J 



{z 2 d Y3 d Yl y 



[(?2 + 0" 3 ]. 



x Zl x Y* 1 Y° 2 Yp [Gi(y, Z)] v . 



(A.9) 



A. 2 PM equation 

The first step consists in recasting the higher-order PDE: 

£i(/2i - Mi) £i(mi - mi + 2) • • • £i(/2i + Mi) C(fii,fjLi;Y,Z) = 0, 
into the following system of equations: 

A (Ai - Mi) Awi = 
A (mi - mi + 2 ) C[ W -i] = 



(A.10) 



(A.ll) 



k A(Ai + Mi)c = q 



[i] 



Afterwards, one expands the C and C[„]'s as in eq. (A.l) , so that the system (A.ll) 
translates into a set of differential recurrence relations for C$ and C^ ai : 



(y 2 d Z3 - y 3 d Z2 ) eg* CT1 = (y dy 2 - y 3 5y 3 + ^) <3£S 



(y 2 9z 3 -y 3 5z 2 ) cf 



[mi-i]o"i 



Y 2 

}'> <>)■., -Y 3 Oy.. 



3 1 2 7 ^[nil'i 

C% i 2 J °[Ml-l]cri + °[Ml]<Tl 



(A.12) 



I {y 2 a Za - y d Z2 ) eg = (y 2 ay 2 - y 3 dy 3 + &±^) ct 1} + cf^ 
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where CL^ = . The solutions C&* to the above system can be written as 

CW(Ai, Mi; Y, ^) = E E fa-a^i Ci(ft + Ml " 2 0; y, Z) , (A. 13) 

where the 's are the expansion coefficients of C [n ] CTl , satisfying 

A(s)C w<ri (x;Y,Z) = 0. (A.14) 

Let us stress once again that, since the C[ n]tri 's may involve singular terms, the expression 
(A. 13) is only a formal solution. Hence, for any (formal) solution C ai to the PM equation, 
there exist (formal) solutions C [n]0 - 1+n to the massless equation such that 

pi n 

C ai (Mi, Mi; ^^) = EE(- 1 fL C w CT1+n (Mi + Ml - 2^; y, Z) . (A. 15) 

n=0£=0 ^ ' 

B Shift solutions 

In this Appendix, we prove that the function (3.52) is a solution to the PM equation (1.8). 
For that, let us first generalize the identities (3.49) to 

C\{x — n) • • • C±(x + n — 2) L\{x + n) Z 2 

= [Z 2 £i(x + n) - (n + l)y 3 ] d(x-n) ■■■ C x {x + n-2), 

L\ [x — n) C\ (x — n + 2) • • • C\ (x + n) Z3 

= [Z 3 £ 1 (x-n) + (n + l)Y 2 ]Ci(x-n + 2) ••• d{x + n). (B.l) 

One can see that the effect of Z 2 (or Z3) is to remove the differential operator L\ with 
the largest (or the smallest) argument. Hence, considering generic powers Z 2 T2 Z^ 3 with 
T 2 + T 3 < Ml) an y solution to the differential equation: 

A(Mi-Mi + 2r 3 ) ••• £ 1 ( / ui + Mi-2r 2 )C T2T3 (y,Z) = 0, (B.2) 

is also a solution to 

A (Mi " Mi) • • • A (Mi + Mi) Z? C T2T3 (y Z) = 0. (B.3) 

Moreover, as any solution to a single L\ equation: 

A (Ml -Mi + 2r 3 + 2n)C T2T3 (y,Z) = 0, (B.4) 

also solves the equation (B.2) for any n = 0, 1, . . . , Ml — T 2 — 73 , one can choose the value 
of n which minimizes |Mi — Mi + 2 13 + 2 n\ . This corresponds to the solution (3.52). 
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C H solutions 



In this appendix we show how to turn a generic function of H into a function of Y and 
Z after integrating by parts all the total-derivative terms present in (3.14). Let us start 
considering the integration by parts of a single Hi : Hi K(Y, Z, H) ~ Hi K(Y, Z, H) . Here 
the H^s are operators defined as 



H, 



i ■ - Yi-i Yi + \ — Ni Zi , 



where 



Ni := Y i+1 d Yl+1 + IS-i dy^ -Yidy t + Zi d Zi + H^H+ph^i 



(C.l) 



(C.2) 



and ~ means equality under the integral sign and modulo TT. The above identity suffices 
to integrate by parts any function of H as 

K(Y,Z,H)= H^H^H^K hlh2h3 (Y,Z) ~ £ H^ H^ H^ K hlh2tl3 (Y, Z) . 

hi,h 2< h 3 h\,h,2,h3 



The operators Hi and Hj commutes when i ^ j , while H i 1 gives 



hi 

£ 

k=0 



[N + h 



ilk 



ZiY(Yi +l Yi 



i-lj 



(C.3) 
(C.4) 



Using the above identity, one can recast eq. (C.3) into a compact form as 

3 



K(Y,Z,H) 



no 



L 



Zid Hi 



Ni+Hid Hi 



K(Y,Z,H) 



H i =Y i - 1 Y i+1 



(C.5) 



In the following, we find a set of functions of H which explicitly give all possible leading 
terms. 



One massless two massive case We aim at finding the functions: 

K(Y 2 ,Y 3 ,H 2 ,H 3 ,Z 1 ), 



(C.6) 



that, after integration by parts, have leading terms involving . Starting from the 

identity (C.5) and requiring that the first n + 1 leading terms cancel, one gets 



Z rj ya 2 ya 3 jjh 2 jjh 3 x 

x(Y 2 H 2 -Y 3 H 3 )J2 (?) 

l — n \ / 



(x - n) k [x + n] n _ k (Y 2 H 2 ) k (-Y 3 H 3 



■ S\ n + 1 



[x + n] 2n+1 Z? Y x h * +h * Y 2 ^ +hs Y 



[x + n] 2n+1 Z? Hp H 3 h3 Yp Y^ [G^Y, Z)\ 



in+1 



TL -yh^+hz va2+h 3 \r(T3+h2 



[Gi(y,z)] n+1 + o(^ 



in+2\ 



(C.7) 



Here, [a] n := a(a — 1) ■ ■ ■ (a — n + 1) and (a) n := a(a + 1) • • • (a + n — 1) are the descending 
and ascending Pochhammer symbols respectively, while x = (u 3 — fJ, 2 )/2 + a 2 — 03 . Let 
us notice that the results come with factors which vanish when x £ {— n, —n + 1, . . . , n} . 
Hence, one has to normalize the corresponding leading terms in order to have non- vanishing 
couplings also for these values of x . 
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Two massless — one massive case In this case, we look for the functions: 



K(Y 3 ,H 1} H 2 ,H 3 ), (C.8) 
which give rise to leading terms proportional to G n+1 , and we get 
Yp H^ 1 H? 2 H% 3 x 

x {H X H 2 - Y 3 2 H 3 ) (x-2n- l) k [x - l] n _ k (H, H 2 ) k {-Y 3 2 H 3 ) n ~ k 

~ (|) n+1 [ x - i] 2n+1 Y 1 ha+ha Y 2 hs+hl Y 3 a:i+hl+h2+n+1 [G(Y, Z)] n+l + 0{5 n+2 ) , (C.9) 
with x = /X3/2 — 03 . 

One PM two massive case In the PM case, a generic Q coupling is of the form: 

mi 

J2y 1 1 k^(y 2 ,y 3 ,h 2 ,h 3 ,z 1 ). (CIO) 

01=0 

On the other hand, the leading terms of the corresponding PI solutions, containing a given 
number n of Z 2 or Z 3 (which can be obtained by using the identity (A. 9)), are 

(Z 2 d Ys - Z 3 dy 2 ) k G?- k Yp +k Y° s+k Z T X \ k = 0, . . . , min{n, Ml } , (Oil) 

where the additional inequality 

a 2 + a 3 + n> a Xl (C.12) 

is enforced. In general, it is rather involved to find a combination of H giving rise to each 
of the leading terms (Oil) since in this case the solutions will be of the form 

z n yn-h yp-hz ya 3 -h+h 3 gh-hs jjh 3 v (k) (y 2 ^ 2 y 3 ^ ^ jr y^ ^ c>13 j 

~ Zl x {Z 2 dv 3 - Z 3 dv 2 ) k G[ l - k Y* 1 Y° 2+k Y 3 a:i+k + 0(5 n+1 ) • (C.14) 
Here, we have introduced the combination: 

V^(Y 2 H 2 ,Y 3 H 3 ,Y 1 Y 2 Y 3 ):= ]T (Y 2 H 2 ) a (Y 3 H 3 ) b (Y 1 Y 2 Y 3 ) c , (C.15) 

a-\-b+c=n 

(k) 

where the coefficients p abc have to satisfy a system of linear equations. Let us notice that 
the inequality (C.12) is automatically satisfied due to the form of H . For a given order 
n in Z 2 and Z 3 , one can count n + 1 different polynomials in Y 2 Z 2 and Y 3 Z 3 . Hence, 
the total number of polynomials with at most n powers of Y 2 Z 2 and Y 3 Z 3 is ( n+1 K n + 2 ) ^ 
matching the number of different coefficients pj¥, entering the combination (C.15). This is 
equivalent to saying that, allowing enough powers of Yy, any leading term in (Oil) can be 
reproduced after integration by parts by a suitable choice of the coefficients in eq. (C.15). 
In the PM case, one has to impose also the bounds c < fii in eq. (C.15) and k < \x\ for the 
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leading terms of eq. (C.ll). Hence, the counting of available leading terms with at most n 
powers of Y 2 Z 2 and Y 3 Z 3 is 



(number of leading terms) 



(Ml+1) 2 (m+2) + (Mi + l)(n- m ) 
|0*i + l)(2n + 2-Mi), 



(C.16) 



while the number of free coefficients in (C.15) satisfying c < is 



(number of pjj>) = ^{m-i + l) = \ ( W + l)(2n + 2 - Ma ) . (C.17) 



mi 



Therefore, as the number of solutions matches, using completeness of the PI solutions 
when (1.10) is not satisfied, one can conclude that the Q solutions span the same space. 

In order to clarify the latter discussion, we provide the explicit example of n = 2 . The 
combination proportional to Gf is given by 



-2( 
-Ay 



[[x + y + l] 3 (x - y + 1) O2H2) 

2 



2 (x + y) 2 [x - y} 2 (Y 2 H 2 ) (Y 3 H 3 ) + (x + y - 1) (x - y - 1) 3 (Y 3 H 3 



(x + y) 2 {x - y) (Y 2 H 2 ) -( x + y )[ x - y} 2 (Y 3 H 3 )\ (Y 1 Y 2 Y 3 ) 
+ 2y (2y - 1) (x + y + 1) (x - y - 1) (Y 1 Y 2 Y 3 ) 2 } , 
while the one proportional to G± (Y 2 Z 2 — Y 3 Z 3 ) is 



(C.18) 



n vm4 y U2-/13 v J 3 _ W3 ijh—hz 



Z? Y 1 



3 YS 



m 



Hi 



x { [x + y + 1] 3 (x - y + 1) (Y 2 H 2 f - (x + y - 1) (s - y - 1) 3 (I3P3) 2 

- 2y [(x + y - 1)3 (x - y) (Y 2 £2) - (x + y) (x - y - 1) 3 (I3 ^3)] Pi Y 2 Y 3 ) 

+ 2 (2y - 1) x (x + y + 1) (x - y - 1) (F x Y 2 Y 3 ) 2 } • (C.19) 

Finally, the combination proportional to (Y 2 Z 2 — Y 3 Z 3 ) 2 reads 

z ti y^i-h y 2 oa-h3 y 3 oa-h+h3 /jh-ha £fts I [a; + 2/ + 1] 3 (a; — j/ + 1) (Y 2 # 2 ) 2 
+ 2 (x + y) 2 [x - y] 2 (Y 2 tf 2 ) (Y 3 tf 3 ) + (x + y - 1) (x - y - 1) 3 (Y 3 H 3 f 
- 4y [(x - 1) (x + y) 2 (x - y) (Y 2 H 2 ) + (x + 1) (x + y) [x - y] 2 (Y 3 #3)] (Y *2 Y 3 ) 

+ 2 (2 x 2 - y) (x - y - 1) (x + y + 1) (F x F 2 y 3 ) 2 } . (C.20) 
Here, x and y are defined as 

(C.21) 



^ + a 3 -a 2 



y = l r-<7i 



As one can see from the above example, if no bound on the powers of Y\ is imposed in 
(C.10), then all possible leading terms can be reproduced. However, for given values of y , 
only a subset of the above solutions satisfy c < \x\. For instance, when [i\ = one can see 
that only the leading terms proportional to G™ are leftover. 
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D Flat-space limit 



In this appendix, we present another argument for the equivalence between the Q solutions 
and PI solutions. It relies on the limiting process where both L and \x go to infinity in 
such a way that their ratio stays finite: 

lim M = M. (D.l) 

Let us notice that, as the supplemental Q and P2 solutions have to satisfy the conditions 
(1.10), they cannot be taken into account in this limit so that effectively we are restricting 
the attention to the Q and PI solutions only. In the (D.l) limit, the PDE (3.29) reduces 
to 

i 1/11+1 

Y 2 d Z3 - Y 3 d Z2 + | (M 2 - M 3 ) d Yl \ C(Y, Z) = 0, (D.2) 

which, for [i\ = , corresponds to the usual flat Noether procedure equation for one mass- 
less and two massive fields. However, in general, it does not have a direct physical inter- 
pretation. The PDE (D.2) can be easily solved as 

Mi 

C = E Y l 1K<Tl 02,^3, Z h H 2 ,H 3 ), (D.3) 

<ri=0 

where the H^s are given by 

H t = Y i+1 + I (Mi - M i+l - Zi , 

= Y l+1 -\d x - [d Xi - d Xi+1 - d Xi _,) Z t . (D.4) 

The (A)dS counterpart of (D.4) can be obtained by adding proper total-derivative terms 
to the Ifcti's, ending up with Hi (3.14). In this way, keeping the number of independent 
solutions unchanged, one recovers the following (A)dS vertices: 

mi 

C = Y. Y i aiR(71( - Y ^ Y 3> Z ^H2,H 3 ), (D.5) 

cri=0 

in agreement with eq. (3.56). Notice that, in the limit, the combinations (C.7) and (C.9), 
giving rise to G™- and G n -type leading terms, become proportional to 

(Y 2 H 2 -Y 3 H 3 ) n , {^Hz-Y^Har, (D.6) 

respectively, that are nothing but the corresponding flat combinations. 
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